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Abstract 



For interactions involving two or more nucleons it is useful to divide pions 
into three classes: potential, radiation, and soft. The momentum threshold for 
the production of radiation pions is Qr = \/WNrn^. We show that radiation 
pions can be included systematically with a power counting in Q^. The leading 
order radiation pion graphs which contribute to NN scattering are evaluated 
in the PDS and OS renormalization schemes and are found to give a small 
contribution. The power counting for soft pion contributions is also discussed. 
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Effective field tlieory is a useful method for describing two-nucleon systems |T[ . Recently, 
Kaplan, Savage, and Wise (KSW) introduced a power counting which accounts for the 
effect of large scattering lengths. A similar power counting is discussed in Ref . 0] . According 
to the KSW power counting, the leading order calculation involves only a dimension 6, four- 
nucleon operator which is treated nonperturbatively. Higher derivative operators and pion 
exchange are treated perturbatively. In evaluating diagrams with pions, three types of 
contributions can be identified: potential, radiation, and soft. Since these pion effects differ 
in size, they each have a different power counting. This distinction arises because there are 
several scales associated with two nucleon systems. In this respect the theory is similar to 
NRQCD and NRQED §. 

In NRQCD there are three mass scales associated with non-relativistic systems con- 
taining two heavy quarks: the heavy quark mass M, momenta ~ Mf , and kinetic energy 
~ Mf^, where v is the relative velocity. QCD effects at the scale M are integrated out 
and appear as local operators in NRQCD. The remaining low energy contributions can be 
divided into potential, radiation (sometimes referred to as ultra-soft), and soft pieces p|-p 



Potential gluons have energy of order Mv"^ and momentum of order Mf , radiation gluons 
have energy and momentum of order Mf^, and soft gluons have energy and momentum of 
order Mv. The power counting for radiation gluons requires the use of a multipole expan- 
sion at a quark-gluon vertex The v power counting of potential and radiation gluons 
can be implemented in the effective Lagrangian by introducing separate gluon fields and 



rescaling the coordinates and fields by powers of v . In Ref. |]10l the separation of scales 
was achieved on a diagram by diagram basis using a threshold expansion. The potential, 
radiation, and soft regimes were shown to correctly reproduce the low energy behavior of 
relativistic diagrams in a scalar field theory. In Ref. it was pointed out that these effects 
may be reproduced by an effective Lagrangian in which separate fields are also introduced for 
the soft regime. Note that soft contributions come from a larger energy scale than potential 
and radiation effects. The heavy quark system does not have enough energy to radiate a 
soft gluon, so they only appear in loops. In Ref. it was shown that soft contributions 
to scattering do not appear until graphs with two or more gluons are considered. 

In the nucleon theory there is another scale because the pions are massive. For the pur- 
pose of power counting it is still useful to classify pion contributions as potential, radiation, 
or soft. For a pion with energy go and momentum g, a potential pion has go ~ g^/M where 
M is the nucleon mass, while a radiation or soft pion has go ~ g > fn-w- In non-relativistic 
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theory, integrals over loop energy are performed via contour integration. Potential pions 
come from contributions from the residue of a nucleon pole and give the dominant contribu- 
tion to pion exchange between two nucleons. For these pions, the energy dependent part of 
the pion propagator is treated as a perturbation because the loop energy, go ~ /M <^ q. 
The residues of pion propagator poles give radiation or soft pion contributions. The power 
counting for soft and radiation pions differs. For instance, the coupling of radiation pions to 
nucleons involves a spatial multipole expansion, while the coupling to soft pions does not. 

In section A, a power counting for graphs with radiation pions is given. To illustrate 
this power counting, we compute the leading radiation pion contribution to S-wave nucleon- 
nucleon scattering. The power counting for soft pion contributions is discussed in section B, 
and an illustrative example is given. We emphasize the importance of not double counting 
when adding radiation and soft pion contributions. 



A. Radiation Pions 

In chiral perturbation theory the expansion is in powers of momenta and the pion mass 
m^. For power counting potential pions it is convenient to take the nucleon momentum p = 
Mv ~ so f = uIt^/M ~ 0.15. The situation is different for radiation pions. There 

is a new scale associated with the threshold for pion production, which occurs at energy 
E = m^r in the center of mass frame. This corresponds to a nucleon momentum p = Qr, 
where Qr = ^/MuIt^ = 360 MeV. Because the radiation pion fields cannot appear as on-shell 
degrees of freedom below the threshold E = 171^^, one expects that the radiation pion can 
be integrated out for p <ti Qr- (Potential pions should be included for p > 171^,/ 2.) Another 
way to see that radiation pions require p ~ Qr is to note that in order to simultaneously 
satisfy fcg = fc^ + and fco ~ ^ ~ Mv"^ requires v ~ ^m-j^/M ~ 0.38. 

The full theory with pions has operators in the Lagrangian with powers of which give 
all the dependence. If the radiation pions are integrated out, then the chiral expansion 
is no longer manifest because there will be m^r dependence hidden in the coefficients of 
operators in the Lagrangian. One is still justified in considering the same Lagrangian, but 
predictive power is lost since it is no longer clear that chiral symmetry relates operators with 
a different number of pion fields. Also, the dependence induced by the radiation pions 
may affect the power counting of operators. For example, as shown in Ref. integrating 
out the pion in the one-nucleon sector induces a nucleon electric polarizability Oe cc 1/m-,^. 
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Alternatively, one can keep chiral symmetry manifest by working with coefficients in the full 
theory and including radiation pion graphs. This is the approach we will adopt. 

The presence of the scale Qr modifies the power counting of the theory with radiation 
pions. In the KSW power counting, one begins by taking external momenta p ~ ~ Q. 
The theory is organized as an expansion in powers oi Q. To estimate the size of a graph, 
loop 3-momenta are taken to be of order Q. However, potential loops within graphs with 
radiation pions can actually be dominated by three momenta of order Q,.. To see how this 
comes about, consider as an example the graph shown in Fig. Ic. Let q be the momentum 
running through the pion propagator, and let k be the loop momentum running through a 
nucleon bubble inside the radiation pion loop. The poles from the pion propagator are 

(ll - - "^l + (go - y^g 2 + + ie) (go + \J q"^ + ml - ie) 

so the radiation pion has |go| > m^. This energy also goes into the nucleon bubbles. The 
k integrand is largest when the nucleons are close to their mass shell. But since the energy 
going into the loop is ~ uIt^, this occurs when fc^/M ~ m^r, i.e., k ~ Q^. We will begin by 
considering the contribution of radiation pions to elastic nucleon scattering at the threshold, 
E = tjIt^. At this energy, external and potential loop momenta are of the same size and power 
counting is easiest. Because p ~ it is obvious that we want to count powers of Qr rather 
than Q. 

Before discussing the power counting, recall the Lagrangian with pions and nucleons P| : 

= ^Tr (9^S + 4^ Tr(m,S + m,St) + ^ N^a^.^ - i'^d^ON 

o 4 Z 

+ N^{iD, + Cl,'\N^Pt^Ny{N^Pt^N) (2) 

Here qa = 1-25 is the nucleon axial- vector coupling, E = is the exponential of pion fields, 
/ = 131 MeV is the pion decay constant, = ^{(.rrig^ + ^^"^9^^); where rrig = diag(m„, m^) 
is the quark mass matrix, and = w{mu + m^). The matrices P/'*'' project onto states of 
definite spin and isospin, and the superscript s denotes the partial wave amplitude mediated 
by the operator. This paper will be concerned only with S-wave scattering, so s = ^5*0, '^S'l. 

The Co operator mediates S-wave nucleon transitions. The D2 operator is important 
because it will be necessary to introduce counterterms proportional to to regulate UV 



+ 



(Ar^P/^)Ar)t(iV^p/^) V^AT) + h.c. 
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divergences appearing in the graphs evaluated below. The parameters appearing in Eq. (0) 
are bare parameters which require renormalization. For systems with two or more nucleons, 
it is necessary to introduce finite subtractions in order to obtain manifest power counting. 
Two such renormalization schemes are Power Divergence Subtraction (PDS) 0,^] and a 
momentum subtraction scheme (OS) These renormalization schemes are discussed 

extensively in Ref. WM. The renormalized coefficient CoifiR) ~ Arc / {M ^r) where is the 
renormalization point. A loop with two nucleon propagators gives ~ Mp/ (47r), so for fiR ^ p 
the Cq{^r) bubble graphs should be summed. For a scattering length a, this summation 
includes powers of {ap) to all orders p|, as desired since a is large. The coefficients C2{fiR) 
and D2{^ir) scale as ~ so for yU/j ~ p ~ m-,, they are treated perturbatively. The 

ellipsis in Eq. (j^) denotes higher order terms including contact interactions which for NN 
scattering begin to contribute at next-to-next-to- leading order (NNLO). These will not be 
considered here. 

Next, we introduce the power counting at the scale Q,- A scheme with manifest power 
counting will be used, so that Cq^hr) ~ l/(MyU/j), C2(/U_r) ~ l/(MAyu|j), etc., where A is 
the range of the theory. We will take p ~ /i_R ~ Qr- A radiation loop has qo ^ q ^ so 
d^q ~ Q^/M"^, where q is the momentum running through the pion propagator. A radiation 
pion propagator gives a M'^/Q^, while the derivative associated with a pion-nucleon vertex 
gives Ql/M. A nucleon propagator gives a M/Ql. External energies and momenta are 
kept in the nucleon propagator since E ~ p^ /M ~ Ql/M. Furthermore, it is appropriate 
to use a multipole expansion for radiation pion-nucleon vertices which is similar to the 
treatment of radiation gluons in NRQCD 0. Therefore, radiation pions will not transfer 
three-momenta to a nucleon. This is usually equivalent to expanding in powers of a loop 
momentum divided by M before doing the loop integral. A potential loop will typically have 
running through it either an external or radiation loop energy ~ Ql/M. Therefore, in these 
loops the loop energy /cq ~ Q^JM^ while the loop three momentum k ~ Q^., so dl^k ~ Ql/M. 
It is not inconsistent for k while q ~ Ql/M, since three momenta are not conserved 

at the nucleon-radiation pion vertices. At the scale Qr potential pion propagators may 
still be treated in the same way, V(^o ~ k"^ ~ fn'i) = —i/ik"^ + fn^) + 0{kl/k^), which has 
an expansion in Ql/M"^. We will see through explicit examples that this power counting 
correctly estimates the size of radiation pion graphs. 

Note that only the potential loop measure gives an odd power of Q,., so without potential 
loops the power counting reduces to power counting in powers of m^r. The power counting 
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FIG. 1. Leading order radiation pion graphs for A^A^ scattering. The sohd hnes are nucleons, 
the dashed hnes are pions and 6M, 5Z are the mass and field renormahzation counterterms. The 
fihed dot denotes the Co{fj,ji) bubble chain. There is a further field renormahzation contribution 
that is calculated in text, but not shown. 

here therefore correctly reproduces the usual chiral power counting used in the one nucleon 
sector |1T5| . 

Graphs with one radiation pion and additional higher order contact interactions or po- 
tential pions are suppressed by factors of Qr/A relative to graphs with a single radiation 
pion and Cq vertices. The Qr expansion is a chiral expansion about = 0, so there is 
a limit of QCD where it is justified. The scale A is unknown. One possible estimate is 
Ann = 87r/^/(M(yf^) = 300 MeV since a graph with m + 1 potential pions is suppressed by 
p/300MeV relative to a graph with m potential pions. However, this order of magnitude 
estimate only takes into consideration a partial subset of the graphs of the theory. As argued 
in Ref. it is possible that the range is of order the scale of short range interactions that 
are integrated out, implying A ~ 500 MeV. In fact, the accuracy of NLO computations of 
nucleon-nucleon phase shifts is in agreement with this physically motivated estimate of the 
range. We will assume that an expansion in Qr/A is valid. This hypothesis will be tested 
further by seeing how well the effective theory makes predictions at p ~ 300 MeV. For 
example, processes with external pions could be considered. If the Qr/A expansion is not 
convergent, then application of the theory is restricted to p < Qr- 

The radiation pion graphs that give the leading order contribution to nucleon-nucleon 
scattering are shown in Fig. |I|. The filled dot denotes the leading order interaction between 
nucleons, a Cq{hr) bubble sum. We illustrate the power counting with an example, the graph 
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in Fig. Id. For the moment, replace the Co bubble sums with single Co vertices. Each Co 
gives a factor of 1/MQr and each nucleon line gives a factor M/Q^. The derivatives from 
the pion couplings combine with the radiation pion propagator to give a factor of unity. 
The radiation loop gives Q^/M^, while the nucleon bubble loop gives Ql/M. There is also 
a factor oil/ p from pion exchange, and two factors of l/47r from the radiation loop giving 
a 1/A^. (A^ ~ 1 GeV is the chiral symmetry breaking scale.) Combining all factors, we find 
that this graph scales like Ql/{M^K^. This graph is suppressed relative to the leading order 
amplitude, A^~^\ by a factor of (5^/(M^A^) = m^/A^. Note that Co bubbles are summed 
on external nucleon lines as well as in the interior of the radiation loop, and each graph in 
the sum has the same size. It is straightforward to verify that all graphs in Fig. 1 scale the 
same way. For external bubble sums we can simply use the vertex iA^^^'^ where A^~^'^ is the 
leading order S-wave amplitude, 

A(-^) = --^, (3) 

and the pole 7 = An/MCoifiR) + fJ'R ~ l/ct- Graphs with two radiation pions are suppressed 
by at least Qf/{M'^A^) = m^/A^ and will not be considered. 

The first graphs we consider are those in Fig. |l|a,b. These graphs have contributions from 
potential and radiation pions, and it may not be obvious that a clean separation occurs. Here 
the energy integrals will be evaluated without any approximations, after which the graphs 
split into radiation and potential parts. The graph in Fig. [^a gives: 

^ , (-1) 9l_ f i i ... 

(Throughout this paper we will include a factor of (/i/2)^~^ in the loop measures.) Per- 
forming the go integral gives a term from the residue of the nucleon pole and a term from 
the pion pole, 

,2 ^ Jn„ T\/r 



(5) 



2P J (27r)" (g - p)^ - ME ^2 + < - [f - 

/ S^^L= , ' ..... , ' ... , (6) 



where n = D — 1. Eq. is the potential pion contribution. Expanding in — isULP 



2M 



L 2M 



gives the result in Ref. 0J3]. The subleading terms in this expansion are suppressed 
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byQ rn^/M"^. Eq. (^) is the radiation pion contribution. With |g| < M, we may take 
[q — •p)'^ /M /M in the last two propagators, which is the same approximation that is 
made by performing the multipole expansion. Finally, we use the equations of motion to set 
E —p^/M = 0. It is important to note that we have not neglected E relative to For 
n = 3 — 2e, Eq. (Bl) becomes 



1 1 , /'mi\ 
- + 7T-ln(^ 
e 3 ^ fi'^ r 



(7) 



(27r)" (g 2 + m2 )3/2 

where /I^ = fi'^ne'^'^ . Note that this integral is finite in three dimensions {n = 2). 

The next graph we consider is shown in Fig. Pd. We have chosen to route loop momenta 
so that q runs through the pion and ±k and ±{k + q) run through the nucleon lines. The 
momentum k is potential, while q can be potential or radiation. Doing the ko contour 
integral and combining the two terms gives: 
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Doing the go integral gives two terms, but the radiation and potential contributions are still 
mixed. Combining these gives 



,[^{-1)12^ 



2 

I 

2/2 



d'^k d^q q- 
(27r)2" 



1 



1 



1 



q^ + ml E-^ E-^ E 



2M 



fc2 



(9) 



fx? + < 



which can be split into potential and radiation parts 



d'^k 



d^'q 



q- 



1 



(27r)«/ (27r)"y^^^^ 



2 p _ fc£ 

^ M 



Ve- 



M 



+ 



g 2 + ml 
1 



(2fc-g+92) 
2M 



(10) 



E-M- 



(2fc-g+q2) 
2M 



g 2 + 



The first term in Eq. (0) is the two-loop potential pion graph evaluated in Ref. [Q. The 
factors of {2k ■ q + g 2) / (2M) appearing in the denominators can be dropped because the 



^ Note that rn^/q^ ~ itLt^/M, but we have kept the m?^ term in the potential pion propagator 
in Eq. (^). We could consider expanding in m^r/g using the asymptotic expansion techniques 
discussed in section B, but for p ~ m^r these terms would have to be resummed. Unlike the 
soft and radiation contributions, there is no issue of double counting for potential pions, so 
for simplicity we will simply keep the th?^ in the propagator. 
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loop integral is dominated hj k, q <^ M and therefore {2k ■ q + g*^) / (2M) <^ yq"^ + rn^. For 
the second term, which is the radiation pion contribution, this is equivalent to the multipole 



expansion. Momenta k ~ y/Mm^ and q ~ dominate the integrals in the second term. 
In Ref. 1 10], the potential and radiation parts of the graph in Fig. |l]b were evaluated in the 
limit = 0, and shown to correctly make up the corresponding part of the fully relativistic 
calculation. The calculation in Ref. |jTO| agrees with Eq. ( |IOD for = 0. Note that the 
radiation part would not agree if we assumed k ~ and used static nucleon propagators 
in the radiation loop0. For n = 3 — 2e the radiation part of Eq. ([T0|) is 



b) = t[A^-^^] 



l)l2 9a 



d^k 



2P J {2tiY J (27r)"g2 + ^2 ^ _ ^ 



1) 



^ g\Mml /3p 
' (47r/)2 [47r 



i4-2.„2-,„(5)-l„(Z|:) 



M 
2 



M 



g 2 + ml 



(11) 



where 



h{.x) = - 



2 r(f: 



3-^2 



+ 



4A 



^ 3-^2 



, 3 13, ,37, 2 
'-~4' 4' 4J' '-2' 4-^'^ 



:i2) 



For n = 2 the loop integral in Eq. ([11]) is finite, except for the function Ji which has a 
p^/in — 2) pole. In PDS this pole would effect the running of C2{fiR), but we will see 
that contributions proportional to Ji will cancel in the sum of graphs in Fig. |l|. Since 
^(-1) l/(MQr) the results in Eqs. (0, |Tl|) are order Q^/(M^A2) as expected. At one-loop 
the 1/e pole in Eq. (|^ is cancelled by a countertermQ 



^uv.la^^ = _3 c, 



finite 9a 




jE + In (vr) 



(13) 



(4vr/)2 

For higher loops the 1/e poles in Eq. (|^), Eq. (|11|), and 5"^'^'^D2 dressed with Cq bubbles 
cancel. Note that the 0(e) piece of the bubbles give a finite contribution. 



Mp\ 
17/ 



1 + e 



2-21n2-ln 



—p — le 



(14) 



The result of combining Figs. |l]a,b and 5^^'^°'D2 dressed by Cq bubbles is: 



2Furthermore, if static nucleon propagators are used one obtains a linear divergence requir- 



ing a non-analytic counterterm oc |16 



^The bare coefficients are written as C^^--^ = S""^ C + C^""'^" . In OS and PDS additional finite 
subtractions are made so that (7^°!*'= = C{^r) — Yj^^C{^r), see Ref. |ll4 |. 
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a)+b)= 3t 



-1)12 



(47r/)2 Att 



+ i 



r .(-1)12 9AMml ^Mm^ .E^ 



_3 

Next we consider the graphs in Fig. [I|c,d,e. The loop integrals in these graphs vanish 
if the pion pole is not taken so there is no potential pion contribution. As pointed out 
in Ref. emission of the radiation pion in these graphs changes the spin/isospin of the 
nucleon pair. Therefore, if the external nucleons are in a spin-triplet (singlet) state, then the 
coefficients appearing in the internal bubble sum are 

Co'^"Vi?) (cf^'Vi?))- The notation 
Co (Cq) will be used for vertices outside (inside) the radiation pion loop. We begin with 
Fig. |l]c. The contribution from the graph with m nucleon bubbles in the internal bubble 
sum is 

dPq i i —iq"^ 



9l 



2PJ (27r 



X 



go - «e go 



it ql 



mi + ie 



m+l 



f27r 







go + f 



2M 



+ ie ko + ^-^ + te^ 



(16) 



where we used the multipole expansion and then the equations of motion to eliminate E 
and p from the first two propagators. All nucleon propagators have a go pole above the real 
axis, while the pion propagator has one pole above and one below. Therefore, the go contour 
is closed below. The d^k integrals are also easily performed giving 



P 



-C^{^^R)MT{l-n/2) 
(47r)"/2 



-p^ + MJ g 2 + m. 



,2 ^71/2-1 



(2vr) 



(g 2 + )3/2 



(17) 



Note that the size of the loop momenta k in the nucleon bubbles is ~ \^Mm^ even for 



p < y/Mrn^. The /i/j inside the brackets comes from inclusion of the PDS or OS 5"Co(yUi?) 
counterterm graphs for the internal bubble sum. The integral will be dominated by g* ~ 
so the graph will scale as 



(18) 



1 ml f^/Mrn^ 

A2 MllR \ flR 

Since /i/j ~ ^/Mm^, all graphs in the sum are of order Q^/{M^A'^). 

For Figs. |l]c,d,e the sum over bubbles should be done before the radiation loop integral. 
The reason is that an arbitrary term in the bubble sum has a much different dependence 
on the energy flowing through it than the sum itself. This can be seen in the g depen- 
dence in Eq. ([TtI) . If we integrate over q then terms in the sum may diverge whereas the 
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integral of the complete sum is finite. In fact, for n = 3, Eq. (|T^) has divergences of the 
form r(— 1 — m / A)Y [E"^ / m^^) and r(— 1/2 — m/A)EF[E'^ /mD where F is a hypergeometric 
function. These divergences are misleading because for momenta > 1/a we know that the 
correct form of the leading order four point function falls off as 1/p. For this reason the 
summation is performed before introducing counterterms to subtract divergences. (This 
approach is also taken in the analysis of three body interactions in Ref. [|r^). Summing over 
m, Eq. ([T7|) becomes: 



c) = —i 



1 



— p2 -)- MJ g 2 -|- 77^2 



n/2-l 



(19) 



^ M 

where 7' = An / MCq^hr) + /ii? ~ 1/a. As expected the graph scales as Q^. In the limit 
n — >■ 3, /2 is finite and given by 



hix) 



r(-!) 



3-^2 



sJ_ 1 h ii l\ 2^ ^£(1 

1 4' 4' 4i'^2' A^'^ ) 2 r 



'-~4' 4' 4J' '-2' 4-^'^ 



(20) 



I2 is manifestly /iij independent and is also finite as n — > 2. 

Next we consider the graph in Fig. |I]d. Integrals are done in the same manner as that of 
Fig. |I]c. For n = 3 - 2e, Fig. |l]d is 

-4.A(-i)-^^ r(n/2 



d) 



(^-p2 + My/g^ + m2 



2\n/2-l 



2/2 (4.)"/2-i 7 (2vr)« (g- + ^2)3/2 _ ^ ^ ^ ^ 



7' 



-p 
n/2-1 



(4ir/)' 



1 , /m2 



^ (p - ^V) MA(-i) ^2 3/2 ^ E ^ 

^^7^ — A^jpyji) ^<^J- ^^^^ 



Fig. |l|d is finite for n = 2, except for the function Ii. The 1/e pole in Eq. (^) is cancelled 
by a new tree level counterterm i5^"'^^D2Tn^ where 5^"'^^D2 has the same form as Eq. (p!3D 
except with a —12 instead of a —3. 

Evaluation of Fig. |l]g is also similar to Fig. ^c. For n = 3 — 2e we find: 
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-2z 



. (p - ii) 
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1 7 
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20F 



Hi?) 



(22) 
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This graph is finite for n = 2, except for the function /i. A D2 counterterm cancels the 
divergence in the last line, 



9\ 



'jE + In (tt) 



(23) 



and (J^^'i'^Do 



Att (47r/)2 

For two and higher loops the remaining 1/e poles cancel between Eqs. (0,^,|£O|j aiiu ■ 1^2 
dressed with Cq bubbles, so no new counterterms need to be introduced. The 0(e) piece 
of the bubbles again give a finite contribution. Combining Figs. |I]c,d,e, and 5"^'^'^D2 and 
^uv,ie£^^ dressed with Cq bubbles gives 



c)+d) + e) = 2i [A^-^ 



-l)l2 



9l 



6mt 



M(7-772) 



1 . /m! 



An 

(7 - 7')' (Mm^)3/2 



-Inf^ 



3 ""V 



^3/2^5/2^ ^ 
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(24) 



Fig. shows a two loop graph with a nucleon self energy on an internal line. It is 
important to also include graphs with the one-loop wavefunction and mass renormaliza- 
tion counterterms, 5Z, 5M inserted on the internal nucleon line. We will use an on-shell 
renormalization scheme for defining these counterterms, which ensures that the mass, M, 
appearing in all expressions is the physical nucleon mass. The counterterms are: 



5M 



3g\ml 



5Z 



9 g\ml 



IGnp ' 2(47r/)2 \e 

The result from the graphs in Fig.|l|f is then 



-In 



(25) 



E 



(26) 



(47r/)2 4v/^ 

When Eq. ( pB]) is added to Eqs. (|1^,^) the terms proportional to Ji cancel. To implement 
PDS we must consider the value of the graphs in Fig. |I|f using Minimal Subtraction with 
n = 2. For n = 2 + e we have SM = 3g\m'l^/ (IGvr/^ e) and 6Z = 0, which makes the sum of 
graphs in Fig. |l|f finite. Finally, renormalization of the bare nucleon fields in the Lagrangian, 
jybare _ ^J~2iN ^ Z = 1 + 6Z, iuduccs a four- nucleon term 



sc 



_^{s),finite ^2^^) (iV^p/^)iV)t(iV^P/^)iV) . 



(27) 



Since 6Z Q"^ this term is treated perturbatively. A tree level counterterm 



2 -1 



finite 9a 

° (4vr/)2 



--1E + ln(7r) 



(28) 
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is introduced to cancel the 1/e pole. Dressing the operator in Eq. 



- 9i [A'^-^ 



-l)l2 



with Co bubbles gives 
(29) 



Att (47r/)2 

Again, for n = 2 we have SZ = so no new PDS counterterms were added. Note that if 
we had instead used bare nucleon fields then there would be no correction of the form in 
Eq. p7|). However, Eq. (p6|) would be modified because the last graph in Fig. |l|f is no longer 
present. When this is combined with the contribution from the LSZ formula the sum of 
Eq. ( p6D and Eq. (^) is reproduced. 

For PDS, the graphs in Figs.|l|a-f are finite for n = 2 so no new finite subtractions were 
introduced. For n = 3, counterterms are introduced to renormalize the terms with ln(//^) 
in Eqs. (|T5|,^,^D|) (in PDS /i = fifj). In OS only terms analytic in are subtracted 



(including m^ln(/i^)). We find D2{jJiR) 



AD2{fiR) = 6Co{^lR) 



9l 



^ D2{fiR) + AD2{fiR), with 

M(7-y) 



An 



-i + . + ln(4) 

3 v^^iy 



(30) 



(47r/)2 

Here k = 1/3 in PDS and k = in OS, and /iq is an unknown scale. Note that the logarithm 
in Eq. (|30|) gives a contribution to the beta function for D2{fiR) of the form 



Mrad) ^ Mil - Y) 

4^2J2 4^ 



R) 



(31) 



This disagrees with the beta function of Ref. 0, because in that paper the beta function 
was calculated including only the one-loop graphs. 

Adding the contributions in Eqs. (1^,24,2^,2^) gives the total radiation pion contribution 
to the amplitude at order Q^: 



iA 



rad 



^ [A-nff Att 



«: + ln(^) -^[A('' 



-l)l2 



AD2{iir) ml 



R) 



+1 



)12 



M(7-7')]2 g\ / E 



Ati 



(32) 



^ M ^ 

The first term here has the same dependence on the external momentum as an insertion of 
the D2 operator dressed by Co bubbles. Its fiR dependence is cancelled by fiR dependence in 
/\D2{^ir)- Note that due to cancellations between graphs, this term is actually suppressed 
by a factor of '~f/Qr relative to what one expects from the power counting. The second term 
in Eq. (^) has a nontrivial dependence on E and is suppressed by an even smaller factor 
of 'j'^/Qr- These cancellations were not anticipated by the power counting and it would be 
interesting to determine if the suppression by factors of 7 ~ 1/a continues at higher orders. 
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FIG. 2. Examples of order Qf radiation pion graphs for NN scattering. 



If so, this might be a consequence of an additional symmetry of the theory in the limit 
a — s> oo. If not, terms at order may actually give the leading contribution of radiation 
pions to NN scattering. 

If we now consider momenta p ~ Q Qr, v^e should fix /i/j at the threshold, 

fiR = \fMrn^, and expand in E/m.„ giving l2{E/mT,) = —3.94 + 0{E/m.„). Therefore, 
the dominant effect of the graphs that occur at order is indistinguishable from a shift 
in D2{fiR). Integrating out the radiation pions amounts to absorbing their effects into the 
effective D2 in the low energy theory. The result in Eq. ( P^ ) is suppressed relative to the 
NLO contributions in Ref. [§] by a factor of roughly 2(7 — '~f')/M ~ 1/10. Since this is 
smaller than the expansion parameter, Q/A ~ 1/3, it can be neglected at NNLO. 

Recall, that in evaluating the graphs in Fig. |l| a mulitpole expansion is used, which in 
this case is an expansion in v = \Jm.,^/M. If the first correction in the multipole series were 
to multiple a nonzero term of order then this would give an order contribution. 

However, we have checked that for all the graphs in Fig. |l]the first correction in the multipole 
expansion gives a vanishing contribution. 

At order graphs such as those in Fig. ^ will contribute to NN scattering. The graph 
in Fig. 0a includes an insertion of the operator 

C = iG2 [N^P^'^N]^ [N^Pt^a^iid.e ~ i'^d.ON] + h.c. . (33) 

(Note that due to the hermitian conjugate this operator is the same for s = ^Sq and s =^5*1.) 
This graph will be dressed with Cq bubbles inside and outside the radiation pion loop. The 
renormalization group equation for G2 gives G2{hr) ~ l/(M/i|,) ~ \/{MQl). Combining 
this with the remaining factors of Qr we find that Fig. |^a is of order Q^/(M^A^) and is 
therefore suppressed by Qr/M relative to a graph in Fig. |l]. Power counting the graphs 
in Fig. 0b, c gives Q^/(M^AnnA^), giving a factor of Qr/ANN relative to a graph in Fig. [l|. 
This provides an example of how graphs with potential pions seem to restrict the range of 
the effective field theory to Ann ~ 300 MeV. The 300 MeV scale applies only to a subset of 
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graphs and may change once all graphs at this order are included. 



It is important to understand how a correction scales with Q for p ~ m-,^. In Ref. 18 



it is demonstrated that graphs can give a Q^^"^'^ contribution plus terms higher order 
in Q. For the Ql calculations the leading terms scale as Q^^"^ and come from the graphs 
in Fig. ^,d,e. However, it turns out the these Q^/"^ terms cancel. It is not understood at 
present why this cancellation occurs. Graphs at order may give an order Q contribution 
for p ~ m^r and should be included in a complete NNLO calculation of the NN phase shifts. 



B. Soft Pions 

In this section soft pion contributions will be discussed. We will see that there are graphs 
with non- vanishing soft contributions that should be included for p > m^. In soft loops two 
scales appear, and p = Mv. We will see below that it is necessary to take p ^ Qr when 
power counting graphs with soft loops in order to avoid double counting. In other words 
V should have the same value as in the radiation pion calculation. A soft loop has energy 
and momentum go ~ ~ Qr? so d^g ~ Qj. The mass of the soft pion is smaller than its 
momentum, and is treated perturbatively. Nucleon propagators in a soft loop are static 
(like in heavy quark effective theory, see eg. Ref. |T^) since the loop energy is greater than 
the nucleon's kinetic energy [^. Therefore, these propagators count as 1/Qr- This power 
counting is identical to that proposed in Ref. ^ except powers of Qr are counted rather 
than Q. 

Unlike potential pions, both soft and radiation pion pieces come from taking the pole 
in a pion propagator. Therefore, we must be careful not to double count when adding 
these contributions. This is accomplished by taking p ^ Qr when evaluating both soft 
and radiation pion graphs. This ensures that the soft and radiation modes have different 
momenta (~ Qr and ~ m^r respectively). Integrals involving the scales Qr and can be 



separated using the method of asymptotic expansions and dimensional regularization |TO| , pO 
Consider splitting a loop integral into two regimes by introducing a momentum factorization 
scale L such that < L < Qr- After the pion pole is taken in an energy integral over qq, 
the remaining integral is of the form 

cTq = J cTq (radiation) + cTq (soft) , (34) 

which is obviously independent of L. In Eq. (R¥) the power counting dictates that expansions 
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a) . . b) , . c) . . d) 



FIG. 3. Examples of one-loop graphs which have soft pion contributions. Graphs a)-d) also 
have a radiation pion contribution, while in addition graph a) has a potential pion contribution. 

in rn^/Ql should be made so that each integral becomes a sum of integrals involving only one 
scale (m^r for radiation and Qr for soft). In dimensional regularization power divergences 
vanish, while logarithmic divergences show up as 1/e poles. Therefore, after expanding 
we can take L — oo in the radiation integral and L — in the soft integral. Taking the 
L oo and L ^ limits may introduce ultraviolet divergences for the radiation integral and 
infrared divergences for the soft integral. When we add the radiation and soft contributions 
any superfluous 1/e poles will cancel. This will be illustrated with an explicit example below. 
The asymptotic expansion procedure has been rigorously proven for Feynman graphs with 
large external Euclidean momenta and large masses It has also been shown to work for 
the non-relativistic threshold expansion of one and two-loop graphs [|r^ . 

Notice that it is crucial to expand the soft pion propagator in powers of m^/Ql, because 
otherwise the radiation pion contribution may be double counted. As an example, consider 
the graph in Fig. ^a. Taking p ^ Qr implies Mv"^ ~ rrij^. For the radiation pion contribution 
go ~ ~ Mf^ ~ m^, so we keep the in the denominator of Eq. (^. When computing 
the soft contribution, we assume qo ^ q ^ Qr ^ and must expand the denominator 
in powers of rn^/Ql. The g integration is now scaleless so the soft contribution to Fig. ||a 
vanishes in dimensional regularization. If we did not expand in m.,^/Qr when evaluating 
the soft contribution, we would have double counted the radiation contribution. The same 
argument can be applied to all the diagrams in Fig. ^ In each case the soft contribution 
vanishes. 

Examples of graphs which have a non-vanishing soft contribution are shown in Fig. |^. 
These diagrams were calculated in Ref. |^ (although the S-wave channels were not analyzed 
there). In the KSW power counting they must be dressed on the outside with Cq bubbles. (If 
Fig. ^a or ^3 are dressed on the inside with Cq bubbles then the soft contribution vanishes.) 
To see how these graphs obtain contributions from the soft and radiation regimes consider 
Fig. I^a. Unlike the massless case |jTO|, this graph has a radiation contribution. In the ^5*0 
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channel the loop integral for Fig. ^ is 



2P 



(27r)^| + go- M 



X 



2M 



+ie f -go-(M+,e ql 



mi + ie 



(-7^2 2--' (35) 
[q — ty — + le 

where t = p' — p and ±p and ±p' are the incoming and outgoing nucleon momenta. Unlike 
the graphs in Fig. ^ we are forced to route an external momentum, t, through a pion prop- 
agator. Taking a nucleon pole in Eq. (|35|) gives the potential pion contribution proportional 
to Mp. Taking the contribution from the pion poles gives soft and radiation contribu- 
tions. Our power counting tells us that the leading order soft contribution will be ~ Q^, 
while the leading order radiation contribution will be ~ Qf/M"^. In Eq. (^5]) the factors of 
E/2 — (g + p)^/(2M) can be dropped. In the soft regime the factors of E/2 — (q + p)'^/ (2M) 
are order Ql/M, and are dropped relative to go ~ Qr leaving static nucleon propagators. 
In the radiation regime E/2 — {q + p)^/(2M) after using the multipole expansion and 
equations of motion. This leaves 

i r qIV r dJ'q q^iq-tf f 1 1 1 



2V2f 



(27r)" _ f)2 I [^-2 + ^2] 3/2 

[q^ + ml]^/Hp-{q~ty 



+ 



1 



q 



It 



(36) 



2V2/V y (27r) 

where n = D — \. The singularity at g ^ = (g — t)^ is cancelled only in the sum of terms 
in the first line of Eq. ([36|). These terms can be calculated separately by introducing an ie 
in this denominator ||10|, giving an average over ±ze as indicated^. The factor of (q — t)'^ 



in the numerator can be removed by partial fractioning. For the soft contribution the scale 
of the loop momentum is set by the external momentum, go ~ |g| ~ |i| ~ Qr- Expanding 
Eq. ( p6D in m^/g^ gives 
rf"g 



\2f 



E 



r(-i/2) 



1927r2 V2/ 



(27r)" ^2q- 1- t ± te) r(-l/2 - m)V{m + 1) V g 

- 18 ml) - - + ln 



9\ 



18 

ml 



1 , fJJ? 







Vg^ 


m 




4 

ml 


(^90 



140^ + 



+-t^-36m^ + 280^ + 



(37) 



^The second line of Eq. ( p6D is more easily split into soft and radiation contributions. If we 
had used the integrand on the first line we would also have to consider (g — t) ^ ~ m'^. 
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where we have kept the first few terms in the expansion. The soft contribution starts at 
order Qf. as expected. The first 1/e pole in Eq. ( pTj ) is an ultraviolet divergence, while 
the second is an infrared divergence. For the radiation contribution go ~ lo] ~ "^tt ^ 1^1- 
Expanding in (2t ■ g) / t ^ gives 



— ? 



2/2 
—i 



1927r2 V2/ 



d'^q 1 

(27r)" [g 2 + ^2]3/2 
1 



E 



2q-t 



9 A 



12 mi 



+ In 



+ 



1 



+ In 



m2 



90^-140^ + 



m 



482 m: 



6 



24 m: + 39 ^ - — ^ + . . . 



(38) 



^2 3 t4 

The radiation contribution starts out as order Qj/M"^ as expected. Note that only powers 
of = Ql/M appear. The 1/e poles in Eq. ( ^81) are ultraviolet divergences. When the 
soft and radiation contributions are added the infrared poles in Eq. (|37|) cancel a subset of 
the ultraviolet poles in Eq. (|38D . Adding Eq. (^) and Eq. ([38|) we find 

-90 /r,2x . r,2 



— z 



1927r2 V2/2 



+ : 



In (90^-140^ 



ml 



18 In 



t2 



72 In 



(A) 



60 



t2 



t4 



+ . . . 



358 mt 



39 ^ + f + 

t2 3 



(39) 



where the remaining ultraviolet 1/e poles are cancelled by counterterms for C2 and D2. 

If we are interested in making predictions for p ~ m^r, then powers of m^/t^ must be 
summed up. Summing the series in Eq. (B^) gives 



3a) 



—I 



90 mt 



- + \n 

e 



58m, 



(128m^ + 16m2t2 _t 
tJt'^ + Ami 



In 



f2 _^ 4^2 _ f 
f2 _^ 4^2 _,_ ^ 



(40) 



Since the coefficients in the series in Eqs. ( p7| , |38D diverge for p ' = p, Eq. (^0]) should be 
used when integrating over dt to obtain the ^5*0 partial wave amplitude (even for p ~ Qr). 
Eq. ( ^OD agrees with the result of evaluating Eq. ( p6D exactly. Although the asymptotic 
expansion is not necessary for evaluating Fig. it allows us to identify the radiation and 
soft contributions to this graph and verify that the power counting for each regime works. 
We also see that adding soft and radiation pions reproduces the correct answer without 
double counting. Recall that power counting with p ^ Qr was necessary to avoid double 
counting for graphs like those in Fig. |l|. 
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For p ~ Qr the diagrams in Fig. ^are order Ql/i^pA"^), and are larger than the order 
graphs with a single radiation pion in Fig. |l|. Decreasing p to p ~ the graphs in Fig. |^ 
give contributions of the form 



where F is a function. For p ~ m.,^ the graphs in Fig. |^ are order which is smaller than 
the graphs in Fig. |l] which include order my^, uIt^, and terms. It is interesting to note 
that the relative importance of these graphs changes with p. The graphs in Fig. ^ dressed 
by Co bubbles give a contribution that is the same size as a four nucleon operator with 6 
derivatives, C^^jiRip^ , and are N^LO in the KSW power counting. 

It would be nice to see the expansion used in evaluating the radiation contribution to 
Fig. 1^ implemented at the level of the Lagrangian. It is not clear to us how to do this at 
the present time. In the radiation regime, the pion whose pole is taken can be thought of 
as a radiation pion. However, the other propagator gives factors of q- tjt^^ etc., which 
look more like insertions of non-local operators than the propagator of a field. Also, since 
in general p ^ p' ■, the couplings for this second propagator change the nucleon momenta 
and therefore do not involve a multipole expansion. Finally, the result in Eq. (^Op does not 
have an expansion in i?/m^. So unlike the radiation pion contribution computed in section 
A, this contribution cannot be integrated out for p < y/Mrn^. For these reasons, the use of 
the term radiation for this contribution differs somewhat from the usage in section A. 

To summarize, we have introduced a power counting in factors of Qr = y/Mrn^ appropri- 
ate for graphs with radiation pions. The order radiation contributions to NN scattering 
were computed and found to be suppressed by inverse powers of the scattering length. Soft 
pion contributions also have a power counting in Q^. For p ^ Qr they are ~ Q^, but are 
higher order than the radiation contributions for p ~ m^r. Higher order corrections are 
suppressed by factors of Qr/A and whether or not this expansion is convergent is an open 
question. If the range of the two-nucleon effective field theory with perturbative pions is 
really 300 MeV, and the suppression by factors of 7 does not persist at higher orders then 
contributions from radiation pions induce an incalculable error of order m^/A^ to the NN 
scattering amplitude in this theory. The validity of the Qr/A expansion can be tested by 
looking at processes at p ~ 300 MeV such as those with external pions. 
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